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Abstract 

We calculate electromagnetic neutral pion production off three-nucleon bound states (^H, ■'He) at threshold 
to leading one-loop order in the framework of chiral nuclear effective field theory. In addition, we analyze 
the dependence of the nuclear S-wave amplitude £0+ on the elementary neutron amplitude E^" which in 
the case of ^He provides a stringent test of the prediction based on chiral perturbation theory. Uncertainties 
from higher order corrections are estimated. 
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1. Introduction 

Threshold neutral pion photo- and electropro- 
duction off the nucleon is one of the finest reac- 
tions to test the chiral QCD dynamics, see for 
a recent review. While elementary proton targets 
are accessible directly in experiment, pion produc- 
tion off neutrons requires the use of nuclear targets 
like the deuteron or three-nucleon bound states like 
■^H (triton) or ^He. For a recent review on reveal- 
ing the neutron structure from electron or photon 
scattering off light nuclei, see |2]. Of particular in- 
terest is to test the counterintuitive chiral perturba- 
tion theory prediction (CHPT) that the elementary 
neutron S-wave multipole E^^^ is larger in mag- 
nitude than the corresponding one of the proton. 
This prediction was already success- 
fully tested in neutral pion photo- |5] and electro- 
production off the deuteron [6]. However, given the 
scarcity and precision of the corresponding data, 
it is mandatory to study also pion production off 
three-nucleon bound states, that can be calculated 
nowadays to high precision based on chiral nuclear 
effective field theory (EFT), that extends CHPT to 
nuclear physics (for a recent review, see ^He 
appears to be a particularly promising target to ex- 
tract the information about the neutron amplitude. 
Its wave function is strongly dominated by the prin- 
cipal "s-state" component which suggests that the 
spin of ^He is largely driven by the one of the neu- 



tron. Consequently, in this letter we calculate thres- 
hold pion photo- and electroproduction based on 
chiral 3N wave functions at next-to-leading order in 
the chiral expansion. Experimentally, neutral pion 
photoproduction off light nuclei has so far only been 
studied at Saclay |8, 9] and at Saskatoon ifloilTlll. 

In general, one has three different topologies 
for pion production off a three-nucleon bound state 
as shown in Fig. [1] While the single-nucleon con- 
tribution (a) features the elementary neutron and 
proton production amplitudes, the nuclear correc- 
tions are given by two-body (b) and three-body (c) 
terms. Based on the power counting developed in 
lfl2ll . at next-to-leading order (NLO), only the topo- 
logies (a) and (b) contribute. Here, we will specifi- 
cally consider threshold photo- and electroproduc- 
tion parameterized in terms of the electric £0+ and 
longitudinal Lq+ S-wave multipoles. In particular, 
we will study the sensitivity of the three-body S- 
wave multipoles to the elementary E^^^ multipole, 

taking the proton amplitude E^^^ from CHPT (as 
this value is consistent with the data 1, 1 3i] and a re- 
cent study based on a chiral unitary approach lll4ll '). 



2. Anatomy of the Calculation 

To analyze the process under consideration, we 
calculate the the nuclear matrix element of the given 
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Figure 1 : Different topologies contributing to pion production 
off the three-nucleon bound state (triangle), (a), (b) and (c) rep- 
resent the single-, two- and three-nucleon contributions, respec- 
tively. Solid, dashed and wiggly lines denote nucleons. pions 
and photons, in order Topology (c) does not contribute to the 
order considered here (NLO). 



transition operator O as: 

{M'j\d\Mj}^ {if,M'jP',^4„\d\ifrMjPiJy) , (1) 

where if/ refers to the three-nucleon wave function 
and ky, q„, P^m and P'^j^ denote the momentum of 
the exchanged (virtual) photon, produced pion and 
the initial and final momentum of the 3N nucleus, 
respectively. The "iN bound state has total nuclear 
angular momentum J - 1/2 with magnetic quan- 
tum numbers My for the initial and M'j for the final 
nuclear state. J can be decomposed in total spin 
S = 1/2, 3/2 and total orbital angular momentum 
L = 0, 1,2. The total isospin is a mixture of two 
components, T - 1/2 and 3/2. While the T - 1/2 
component is large, the small T - 3/2 component 
emerges due to isospin breaking and is neglected 
in our calculation. Here, we consider neutral pion 
production by real or virtual photons of a spin- 1/2 
particle - either the nucleon or the and -^He nu- 
clei. At threshold, the corresponding transition ma- 
trix takes the form 

Ma = 2/ £,)+ (Q,T ■ 5") + 2i Lo+ • 5^) , (2) 

with Q,T = Q - (q • ky)ky and = ■ ky)ky 
the transverse and longitudinal photon polarization 
vectors. The transverse and longitudinal S-wave 
multipoles are denoted by £0+ and Lo+, respective- 
ly. Note that Lo+ contributes only for virtual pho- 
tons. 

As explained before, the matrix element Eq. ([TJ 
receives contributions from one- and two-nucleon 
operators at the order we are working. Consider 
first the single nucleon contribution, given in terms 
of the 1-body transition operator O'^. After some 
algebra, one finds 



T ^ ^^0+ ^ T 



+ I" "F' 



'S-V\ 
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where Fj"^^ 
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TIL 



and 



F^jJ^ denote the 



corresponding form factors of the 3N bound state. 

Ft Q,T/L -^M'jMj = 2 <^yl^'i.T/L ■ ^\\Mj)^ , (4) 

which parametrize the overall normalization of the 
response of the composite system to the excitation 
by photons in spin-isospin space. In the above equa- 
tion, cr, (t*,) denote the spin (isospin) Pauli matri- 
ces corresponding to the nucleon /. Furthermore, z 
refers to the isospin quantization axis. 

Using the 3N wave functions from chiral nu- 
clear EFT at the appropriate order, the pertinent 
matrix elements in Eq. (|4]i can be evaluated. Here, 
we use chiral 3N wave functions obtained from the 
N^LO interaction in the Weinberg power counting 
1 15 Q In order to estimate the error from higher 
order corrections, we use wave functions for five 
different combinations of the cutoff A in the spec- 
tral function regularization of the two-pion exchan- 
ge and the cutoff A used to regularize the Lipp- 
mann-Schwinger equation for the two-body T-ma- 
trix. The wave functions are taken from Ref. iUtI 
Tsl] and the corresponding cutoff combinations in 
units of MeV are (A, A) = (450,500), (600,500), 
(550,600), (450,700), (600,700). All five sets de- 
scribe the binding energies of the ^He and ^H nu- 
clei equally well. 

The one-body contributions to the 3N multi- 
poles are given by 
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+ ^0+ 



(5) 



Here, K^^ is the kinematical factor to account for 
the change in phase space from the IN to the 3N 
system, 



A'lN — 



1.092, 



(6) 



with nifi being the nucleon mass and m^j^ the mass 
of the three-nucleon bound state. 



' The consistency of the Weinberg counting for short-range 
operators and the non-peiturbative renormalization of chiral 
EFT are currently under discussion, see the review for more 
details. A real alternative to the Weinberg approach for practical 
calculations, however, is not available. 



nucleus 


^He 




pS+V 


0.017(13X3) 


1.493(25)(3) 


pS-V 

1 J 


1.480(26)(3) 


0.012(13X3) 


ps+v 


-0.079(14)(8) 


1.487(27)(8) 


pS-V 


1.479(26)(8) 


-0.083(14)(8) 



Table 1 : Numerical results for the form factors F!^j^ . The first 
error is an estimate of the theory error from higher orders in 
chiral EFT while the second error is the statistical error from 
the Monte Carlo integration. 



Figure 2: Leading two-nucleon contributions to the nuclear pion 
production matrix element at threshold. Solid, dashed and wig- 
gly lines denote nucleons, pions and photons, in order 



We evaluate the matrix elements for the one- 
body contribution in Eq. (|4]l numerically with Monte 
Carlo integration using the VEGAS algorithm 1,191. 
The results for the form factors Fj^j"^ are given in 
Table [U The first error represents the theoretical 
uncertainty estimated from the cutoff variation in 
the wave functions. We take the central value de- 
fined by the five different cutoff sets as our predic- 
tion and estimate the theory error from higher-order 
corrections from the spread of the calculated val- 
ues. Strictly speaking, this procedure gives a lower 
bound on the error, but in practice it generates a 
reasonable estimate. 

We stress that we follow the nuclear EFT for- 
mulation of Lepage, in which the whole effective 
potential is iterated to all orders when solving the 
Schrodinger equation for the nuclear states. As dis- 
cussed in Ref. [20], the cutoff should be kept of 
the order of the breakdown scale or below in or- 
der to avoid unatural scaling of the coefficients of 
higher order terms. Indeed, using larger cutoffs can 
lead to a violation of certain low-energy theorems 
as demonstrated in Ref. 112 ill for an exactly solvable 
model. 

The error related to the expansion of the pro- 
duction operator is difficult to estimate given that 
the convergence in the expansion for the single nu- 
cleon S-wave multipoles is known to be slow, see 
Ref. f?] for an extended discussion. We therefore 
give here only a rough estimate of this uncertainty. 
The extractions of the proton S-wave photoproduc- 
tion amplitude based on CHPT using various ap- 
proximations 1.22.1 lead to an uncertainty A/iy^ ^ 
+0.05 X lO-^/M^+, which is about 5%. Similarly, 
we estimate the uncertainty of the neutron S-wave 
threshold amplitude to be the same. Consequently, 
our estimate of the error on the single nucleon am- 
plitude is 5%. 

The statistical error from the evaluation of the 
integrals is typically one order of magnitude smaller 



than the estimated theory error and can be neglected. 

We now switch to the two-nucleon contribu- 
tion. In Coulomb gauge, only the two Feynman 
diagrams shown in Fig.|2]contribute at threshold to 
the order we are working [5|. Their contribution to 
the multipoles can be written as 



72N 
"0+ 

f-2N 
-"0+ 



(7) 



with the prefactor 

A'2N 



M„egAmm 



167r(m3N + M„){2n)iFl 
0.135 fm X \Q-^IM„* . 



(8) 



The numerical value for was obtained using 
gA - ^ -26 for the axial coupling constant, F„ - 93 
MeV for the pion decay constant, and the neutral 
pion mass M„ = 135 MeV. The transverse and lon- 
gitudinal form factors F^^^^ and F^*^*^ corresponding 
to diagrams (a) and (b), respectively, are 

(pi2-/l'2 + M2) 

(9) 

and 

^t]l ^-i.T/L ■ ^M'jMj = 3 {M'j\ (ti ■ T2 - r^r\) 
^ [(#12 - Pn - ky/2) ■ ((?! + (?2)] 

[(Pi2 - Pn - kyl2)^ + Ml][pn - Pi^ + Kl^f 
^[eAj^-{pi2-Pn)]\Mj)^, (10) 

where - (k\ - k2)/2 and p'-^2 - i^'i ~ ^2)/'^ 
the initial and final Jacobi momenta of nucleons 1 
and 2, respectively. The integral for the form fac- 
tors F^j.'!^ contains an integrable singularity which 
can be removed by an appropriate variable trans- 
formation. Then, the form factors can be evalu- 
ated using Monte Carlo integration in the same way 
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nucleus 


^He 


3H 


Ff^-Ff^ [fm-i] 


-29.3(2)(1) 
-22.9(2)(1) 


-29.7(2)(1) 
-23.2(1)(1) 



Table 2: Numerical results for the form factors Fi: 



T/L TIL 

parametrizing two-body contributions in units of fm" ' . The first 
eiTor is an estimate of the theory error from higher orders in 
chiral EFT while the second error is the statistical error from 
the Monte Carlo integration. 



One observes that the multipoles get a large con- 
tribution from the two-body terms. This behavior 
is similar to the deuteron case fs*]. For example, 
in the case of Eq+ for ^He, the proton contribution 
is -0.01, the neutron one is 1.72 while the two- 
body contribution is -3.95 in the canonical units of 
10"^/M;r+. There is, however, still a large sensitiv- 
ity to the single-neutron contribution. 

The corresponding threshold S-wave cross sec- 
tion for pion photoproduction a^) is given by 



as the form factors for the single-nucleon contri- 
bution. Our results for F^^'^ - F^^'^ are given in 
Table |2l The first error is again the theory error 
estimated from the cutoff variation in the chiral in- 
teraction as described above. The second error is 
the statistical error from the Monte Carlo integra- 
tion which is about half the size of the theory error 

3. Results and Discussion 

We are now in the position to evaluate the nu- 
clear S-wave multipoles. They are given as the sum 
of the one- and the two-nucleon contributions given 
in Eqs. (|5]|7]) in the previous section. 



Lo+ - 



2N 



+ L; 



0+ 
2N 



(11) 



Using the values for the one- and two-body form 
factors in Tables [T]and|2]together with the sublead- 
ing chiral perturbation theory results for the single- 
nucleon multipoles at 0{p'^) flQl 



Eq+ 



-1.16 X \Q-^IM„. 
-1-2.13 X \Q-^IM^ 

-2.41 X \0-^IM„. 



(12) 



we obtain for the threshold multipoles on ^He and 
on the values in Table [3] For the IN contribu- 
tion, the first error is the theory error from higher 
orders in chiral EFT estimated from the cutoff vari- 
ation as explained above, while the second error is 
from the 5% uncertainty of the one-nucleon ampli- 
tudes. In the case of the 2N contribution, only the 
theory error is given. The total error is obtained by 
adding the theory error and the uncertainty of the 
one-nucleon amplitudes in quadrature. As noted 
before, the total error is dominated by the uncer- 
tainty in the single nucleon amplitudes. We stress 
that our estimate of the theory error is only a lower 
bound. 



uq - 



\ky\ dcr 
\qn\ dQ. 



= |£o+|- 



(13) 



The longitudinal multipole L()+ contributes only in 
electro-production. The corresponding threshold 
cross section contains an extra term ~ |Lo+l^- From 
here on, we will, however, concentrate on photo- 
production. In Fig.|3] we illustrate the sensitivity of 
Go to the single-neutron multipole E^". The inner 
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Figure 3: Sensitivity of oq for He in units of 10 IM\ to the 
single-neutron multipole E'^" in units of 10 IM„+. The ver- 
tical dashed line gives the CHPT prediction for E'^^ and the 
vertical dotted lines indicate the 5% error in the prediction. The 
inner shaded band indicates the theory error estimated from the 
cutoff variation as described in the text. The outer shaded band 
corresponds to a 10% uncertainty in the 2N contribution. 

shaded band indicates the theory error estimated 
from the cutoff variation as described above. The 
outer shaded band illustrates the effect of a 10% un- 
certainty in the 2N contribution. This corresponds 
to the size of the correction to the 2N contribu- 
tion at the next order in the deuteron case ■ The 
vertical dashed line indicates the CHPT prediction 
E"^^ = 2.13 X IQ-'^IM,^. Changing this value by 
+20% leads to changes in a^) of about +30%. Thus, 
the ^He nucleus appears to be a very promising tar- 
get to test the CHPT prediction for E^". On the 
contrary, neutral pion production on ^H is rather 
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nucleus 


multipole 


IN 


2N 


total 


^He 


£(,+ [10-Vm,.] 
Lo+ [lO-^M,.] 


+ 1.71(4X9) 
-1.89(4)(9) 


-3.95(3) 
-3.09(2) 


-2.24(11) 
-4.98(12) 




£0+ [IQ-^IM,.] 
hn [lO-'/M,.] 


-0.93(3)(5) 
-0.99(4)(5) 


-4.01(3) 
-3.13(1) 


-4.94(7) 
-4.12(7) 



Table 3: Numerical results for the threshold multipoles £0+ and Lo+ on 'He and ^^H. The error estimates are explained in the text. 



insensitive to E^^": a variation of E^" from to 3 
changes qq only by 1%. 

Next we compare our predictions with the avail- 
able data. The consistency of the CHPT predic- 
tion for the single-neutron multipole with the mea- 
sured S-wave threshold amplitude on the deuteron 
from Saclay and Saskatoon is well established, see 
Refs. miifl. The reanalyzed measurement of the 
S-wave amplitude for ■'He at Saclay gives £3 = 
(-3.5 + 0.3) X IQ-^/M^. [8, 9], which is related to 
ao according to 

Here, we have approximated the A = 3 body form 
factor of Argan et al. ^ by the numerically 
dominant form factor Fj^^ for ^He, cf. Tab.[T] This 
results in 

Eq+ = (-2.8 + 0.2) X lO-^/M„. , (15) 

assuming the same sign as for our •'He prediction in 
Table [3] In magnitude, the extracted value is about 
25% above the predicted one. Given the model-de- 
pendence that is inherent to the analysis of Ref. [(91], 
it is obvious that a more precise measurement us- 
ing CW beams and modem detectors is very much 
called for 

4. Summary and Outlook 

In this letter, we have presented a calculation of 
neutral pion production off ^H and ^He at thresh- 
old to leading one-loop order for the production 
operator in the framework of chiral nuclear effec- 
tive field theory. We used the chiral wave func- 
tions of Refs. lll5i[T6ll which are consistent with the 
pion production operator to calculate the S-wave 
3N multipoles £0+ and L()+. To this order, the pro- 
duction operator gets both one- and two-body con- 
tributions. Our calculation shows that the two-body 
contributions are of the same order of magnitude as 



the one-body contributions. A similar behavior was 
observed in the deuteron case |5]. 

The theoretical uncertainty resulting from the 
cutoff variation in the employed wave functions ap- 
pears to be small (of the order of 3%). The dom- 
inant theoretical error at this order stems from the 
threshold pion production amplitude off the proton 
and the neutron, which is estimated to be about 5%. 

We explored the sensitivity of neutral pion pho- 
toproduction on ^^He to the elementary neutron mul- 
tipole /iy^" and found a large sensitivity. This makes 
^He a promising target to test the counterintuitive 
CHPT prediction for E^l" |3, 4]. The cutoff varia- 
tion estimate leads to a very small error for the 2N 
contribution. If the error of this contribution is ar- 
tificially enlarged by a factor of 10, the extraction 
of the neutron multipole is still feasible experimen- 
tally. 

We have shown that our prediction for the ■'He 
S-wave multipole £0+ is roughly consistent with 
the value deduced from the old Saclay measure- 
ment of the threshold cross section [9]. A new 
measurement using modem technology and better 
methods to deal with few-body dynamics is urgent- 
ly called for 

There are many natural extensions of this work. 
They include investigating higher orders, pion pro- 
duction above threshold, the extension to virtual 
photons and pion electroproduction, production of 
charged pions, and considering heavier nuclear tar- 
gets such as "^He. Further work in these directions 
is in progress. 
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